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The formalism for analysing the magnetic field distribution in Pauli limited superconductors 
developed earlier is applied to the field dependence of the vortex lattice static linewidth measured 
in Muon Spin Rotation (/iSR) experiments. In addition to writing analytical formulae for the static 
linewidth for the vortex structure in the limit of independent vortices (i.e. moderate magnetic fields), 
we use Abrikosov's analysis to describe the field variations of the static linewidth at the approach 
of the superconductor to metal transition in the limit where the critical field is determined by Pauli 
depairing. 
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It has been proposed in a Letter by Spehling et al.^that 
anomalous variations of the static linewidth as measured 
by Muon Spin Rotation (/xSR) on the heavy fermion su- 
perconductor CeCoIns is due to coupling between su- 
perconducting and antiferromagnetic order. The static 
linewidth was found to increase with applied field in a 
wide range below the upper critical field Hc2 whereas 
the theoretical description available so faJ^ predicted a 
monotoneous decrease in as- Here it is shown that the 
chief effect responsible for such a behavior is a novel 
mechanism that was analyzed befor^ and observed in 
Small Angle Neutron Scattering (SANS) experimentPI^ 
on the Vortex Lattice (VL) of CeCoIns at low temper- 
ature and high magnetic field applied along the c-axis. 
Prior publication of analytical expressions resulting from 
the Ginzburg-Landau formulatiorP, numerical result^ 
based on the Bogoliubov equations were reported. 

Because of a large electron effective mass (to* ~ 
lOOTOe), the diamagnetic screening supercurrents gener- 
ated by the electron Zeeman spin response under a field 
are important and their experimental signature domi- 
nates the usual charge response. As a result the field 
distribution is modified on a distance ~ from the 
center of each vortesP. The existence of these currents 
were anticipated^ in the context of electrodynamics of 
the FFLO (Fulde-Ferrel-Larkin-Ovchinnikov) state. 

The VL static linewidth is defined as 
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where 7^ = 27r x 135.5342 MHz/T is the muon gyro- 
magnetic ratio, h{r) is the component of the local field 
parallel to the applied field H, overline means averaging 
over an unit cell of the VL, and the induction B c:^ H 
in the limit of a large Ginzburg-Landau (GL) param- 
eter K = A/^ (A and ^ are the two length-scales of 
the Ginzburg-Landau theory). This quantity can be ex- 
pressed through a sum involving all order Fourier com- 
ponents Fmn of the field distribution in the VL, 
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The components Fmn are designated the VL form factors 



(FF)pISl in the context of SANS experiments. The mea- 
surement of the first order FF Fio at T = 50 mK irP^l 
has revealed a similar behavior to the one obtained in 
^SR experimenlfi': Fio increases with field up to 4.7 T to 
eventually decrease at the approach of the (first order) 
superconducting to metal transition. 

The FF calculations'^ were performed at temperatures 
where the superconducting transition is second order or 
weakly first order T > 1 K, i.e. where the GL formulation 
is expected to describe qualitatively the vortex lattice 
field distribution. In the limit a/^ ^ 1, with a = yjcj)^/ B 
the inter-vortex distance in a square VL, the FF can be 
written as a sum of two distinct contributions 
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The first contribution is the usual charge response which 
gives rise to the supercurrents we called orbital. In the 
isolated vortex approximation and in the large-K limit it 
writes 
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and shows a monotoneous decrease with applied field. 
Here ^„ = v^C is a variational parameter that mini- 



mizes the total free energy, g„ 
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for a square vortex lattice, Kn{z) is the n"* order modi- 
fied Bessel function of the second kind, and (t>a = 2.07 x 
10^'' G.cm^ is the flux quantum. 

In contrast, the temperature and field dependent Zee- 
man spin contribution is 
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It decreases with temperature, increases with field, and 
brings the main effect at _B < Hc2 and T on the order 
of Tc = 2.3 K. Here /i = (7/is/2 is the electron magnetic 
moment and ^'•^•'(z) is the derivative of the digamma 
function. The above expression is specific to the symme- 
try of the order parameter and is written here for d-wave 
pairing. 

The static linewidth as expressed in Eqs. (3] 5 1 is 
shown in Fig. [l] (here for simplicity we have neglected 
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the small field-dependence of the characteristic lengths 
of the Ginzburg-Landau theory). 



Eq. ( 10 1 shows explicitly that the vortex lattice contribu- 




FIG. 1. Variation of the /iSR static linewidth (Eqs. (2|5l). 
We used the temperature T = 1.3 K, and the parameters 
g = 2 and = 50A. 



The FF given by Eqs. (2][5) is found in the indepen- 
dent vortex approximation. It does not work near the 
transition to the normal state at which point Us drops to 
zero where the transition is of the second order, or to a 
finite value where the transition is of the first order. In 
the high field limit which is close to the transition line, 
the main source of magnetic fi eld inhomogcneity in the 
vortex lattice is the Zeeman on^^ 
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and iVo is the non-superconducting density of states at 
the Fermi level. The Fourier decomposition of the square 
of the gap magnitude reads 
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Therefore the FFs corresponding to the Bragg peaks with 
indices (m, n) ^ (0, 0) take on the form 
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and the vortex lattice static linewidth becomes simply 

or = ^7pe|A(x,y)P, (10) 

where 



+ CX3 

E 



n— — oo 



1 ~ 0.4247. 



(11) 



tion to the static linewidth vanishes when the transition 
is of the second order but shows a discontinuity when the 
transition is of the first order. In the former case, the gap 
average is known'^ ^ 
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Here a and {3 are respectively the quadratic and 
quartic coefficients of the Ginzburg-Landau free-energy 
(a changes sign at the transition), and Pa = 

2 

|A(x, y)|"'/|A(a;, y)P is the Abrikosov parameter. It is 
jf^ = 1.18 for a square vortex lattice and = 1.16 for 
a triangular lattice. Then it follows 
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which is shown in Fig. [2j 
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FIG. 2. Behavior of the pSR static linewidth close to the 
second-order transition line at temperature T = 1.3 K as ob- 
tained from Abrikosov's analysis in the Pauli limit. 

While experiment was made at T = 20 mK, our results 
are shown at T = 1.3 K where the GL formulation can 
be applied. Our model yields an effect that is observable 
in this regime and it is expected to be enhanced at lower 
temperatures. It would be interesting to obtain experi- 
mental information on the static linewidth variations at 
higher temperatures. 

As a final remark, the authors in [Ref. 1] assumed 
a temperature-independent static linewidth to extract 
the temperature dependence of the relaxation rate dy- 
namic contribution. However it is known from the 
temperature dependence of A [Ref.*^] that as varies as 
l-[27rA(0)/r]i/2exp[-A(0)/T] for /xB < TandT < 
(A(0) is the value of the superconducting gap at T = 0), 
and as (T^ - T)/Tc for < T and - T < T^. For 
higher fields, we obtain from Eqs. (2][5) significant and 
qualitatively similar variations. 
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